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QJ Abstract: 



In explaining such phenomena as the redshift and even the fact that we can see stars 
and planets, the effect of wave packet spreading (WPS) of the photon wave function 
is not taken into account. Probably the main reason is a belief that WPS is not 
important since a considerable WPS would blur the images more than what is seen. 
■ However, WPS is an inevitable consequence of quantum theory and moreover this 

^ effect is also known in classical electrodynamics. So it is not sufficient to just say 

that a considerable WPS is excluded by observations. One should try to estimate the 
importance of WPS and to understand whether our intuition is correct or not. We 
^ ■ explicitly demonstrate that a standard relativistic quantum-mechanical calculation 

shows that spreading in the direction perpendicular to the photon momentum is very 
Cif important and cannot be neglected. Hence the physics of the above phenomena is 

not well understood yet. Possible approaches for solving the problem are discussed. 

PACS: 03.65 Sq, 42.15.-i, 95.30.-k, 98.62.Py 

Keywords: quantum theory, semiclassical approximation, redshift 

1 Introduction 

The problem of explaining the redshift phenomenon has a long history. Different 
competing approaches can be divided into two big sets which we call Theory A and 
Theory B. In Theory A the redshift has been originally explained as a manifestation 
of the Doppler effect but in recent years the cosmological and gravitational redshifts 
have been added to the consideration. In this theory the interaction of photons with 
the interstellar medium is treated as practically not important, i.e. it is assumed that 
with a good accuracy we can treat photons as propagating in the empty space. On 
the contrary, in Theory B, which is often called the tired-light theory, the interaction 
of photons with the interstellar medium is treated as a main reason for the redshift. 
At present the majority of physicists believe that Theory A explains the astronomical 
data better than Theory B. Even some physicists working on Theory B acknowledged 
that any sort of scattering of light would predict more blurring that is seen (see e.g. 
the article "Tired Light" in Wikipedia). 
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A problem arises whether or not the effect of wave packet spreading (WPS) 
of the photon wave function is important for explaining the redshift. One might think 
that this effect is not important since a considerable WPS would also blur the images 
more than what is seen. Moreover, the very fact that we can see stars and even 
planets is an indication that for some reasons WPS is not explicitly manifested in 
observational astronomy. However, WPS is an inevitable consequence of quantum 
theory and moreover this effect is also known in classical electrodynamics. Hence 
it is not sufficient to just say that a considerable WPS is excluded by observations. 
One should try to estimate the importance of WPS and to understand whether our 
intuition is correct or not. 

The observational data give grounds to believe that with a very high accu- 
racy we can describe the propagation of light from distant objects in the framework 
of geometrical optics. In turn, it is well-known that the approximation of geometrical 
optics can be obtained from semiclassical approximation in quantum theory. 

The main goal of this paper is to investigate semiclassical approximation 
for ultrarelativistic particles and especially for photons. Our presentation is selfcon- 
tained, we discuss in detail all the assumptions made in the calculation and argue 
that they are fully in the framework of standard quantum theory. Hence our main 
result can be treated as a standard result of relativistic quantum theory in semiclas- 
sical approximation. This result is derived in Sec. HI As a preparatory step, in Sec. 
E] we describe well-known facts about semiclassical approximation in nonrelativistic 
quantum mechanics and in Sec. [3] those results are generalized for wave packets of 
ultrarelativistic particles in momentum space. In Sec. [5] it is shown that in classical 
electrodynamics the results on the validity of geometrical optics can be obtained by 
reformulating the results on the validity of semiclassical approximation in quantum 
theory. Finally, Sec. Elis a discussion. 

2 Semiclassical approximation in quantum me- 
chanics 

In quantum theory, states of a system are represented by elements of a 
projective Hilbert space. The fact that a Hilbert space H is projective means that if 
%l) & H is a, state then const ip is the same state. The matter is that not the probability 
itself but only relative probabilities of different measurement outcomes have a physical 
meaning. In particular, normalization of states to one is only a matter of convention. 
In the present paper we assume this convention, i.e. we will work with states ip such 
that WipW = 1 where ||...|| is a norm. It is defined such that if (...,...) is a scalar 
product in H then \ \ip\ \ = (ip,ipY^'^. 

In quantum theory every physical quantity is described by a selfadjoint 
operator. Each selfadjoint operator is Hermitian i.e. satisfies the property ('02, Aipi) = 
{Ailj2,ipi) for any states belonging to the domain of A. If A is an operator of some 
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quantity then the mean value of the quantity and its uncertainty in state ijj are given 
hy A — {ip, Alp) and AA — ||(^ — respectively. The condition that a quantity 

corresponding to the operator A is semiclassical in state ip can be defined such that 
\AA\ <C \A\. This implies that the quantity can be semiclassical only if \A\ is rather 
large. In particular, if ^4 = then the quantity cannot be semiclassical. 

Let B be an operator corresponding to another physical quantity and B 
and AB be the mean value and the uncertainty of this quantity, respectively. We 
can write AB = {A,B}/2 + [A,B]/2 where the commutator [A, B] = AB - BA 
is anti-Hermitian and the anticommutator {A, B} = AB + BA is Hermitian. Let 
[A, B] = —iC and C be the mean value of the operator C. 

A question arises whether two physical quantities corresponding to the 
operators A and B can be simultaneously semiclassical in state ■0. Since ||'0i||||'02|| > 
I ('01) '02) I , we have that 

AAAB > ^liij, {{A -A,B-B} + [A, B])ij)\ (1) 

Since (■0, {A — A, B — B}t/j) is real and (■0, [A, B]t/;) is imaginary, we get 

AAAB > ^\C\ (2) 

This condition is known as a general uncertainty relation between two quantities. A 
well-known special case is that if P is the x component of the momentum operator 
and X is the operator of multiplication by x then [P, -'^j = —ih and ApAx > h/2. 
The states where ApAx = h/2 are called coherent ones. They are treated such that 
the momentum and the coordinate are simultaneously semiclassical in a maximal 
possible extent. A well-known example is that if 

i^i^) = -^expl^pox - TT^ix- xof] 
ayvr n 2a^ 

then X = xq, P = po, Ax = a/\/2 and Ap = h/{aV2). 

Consider first a one dimensional motion. In standard textbooks on quan- 
tum mechanics, the presentation starts with a wave function ijj{x) in coordinate space 
since it is implicitly assumed that the meaning of space coordinates is known. Then 
a question arises why P = —ihd/dx should be treated as the momentum operator. 
The explanation is as follows. 

Consider wave functions having the form ■0(a;) = exp{ipox/h)a{x) where 
the amplitude a{x) has a sharp maximum near x = xq E [xi,X2] such that a{x) is 
not small only when x & [xi, ^2]. Then Ax is of the order X2 — Xi and the condition 
that the coordinate is semiclassical is Ax ^ |a;o|. Since —ifidip{x)/dx = pQip{x) — 
ihexp{ipox/h)da{x)/dx, we see that iplx) will be approximately the eigenfunction of 
—ihd/dx with the eigenvalue po if |pofl(2;)| ^ h\da{x)/dx\. Since \da{x)/dx\ is of the 
order of \a{x)/Ax\, we have a condition |poAa;| ^ h. Therefore if the momentum 
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operator is —ihd/dx, the uncertainty of momentum Ap is of the order of h/Ax, 
\Pq\ ^ Ap and this imphes that the momentum is also semiclassicaL At the same 
time, \poAx\/27rh is approximately the number of oscillations which the exponent 
makes on the segment [xi,X2]. Therefore the number of oscillations should be much 
greater than unity. In particular, semiclassical approximation cannot be valid if Ax 
is very small, but on the other hand. Ax cannot be very large since it should be 
much less than xq. Another justification of the fact that —ihd/dx is the momentum 
operator is that in the formal limit h ^ the Schroedinger equation becomes the 
Hamilton- Jacobi equation. 

We conclude that the choice of —ihd/dx as the momentum operator is jus- 
tified from the requirement that in semiclassical approximation this operator becomes 
the classical momentum. However, it is obvious that this requirement does not define 
the operator uniquely: any operator P such that P — P disappears in semiclassical 
limit, also can be called the momentum operator. 

One might say that the choice P = —ihd/dx can also be justified from the 
following considerations. In nonrelativistic quantum mechanics we assume that the 
theory should be invariant under the action of the Galilei group, which is a group of 
transformations of Galilei space-time. The x component of the momentum operator 
should be the generator corresponding to spatial translations along the x axis and 
—ihd/dx is precisely the required operator. In this consideration one assumes that 
the space-time background has a physical meaning while, as discussed in Ref. [I] and 
references therein, this is not the case. 

As noted in Ref. [Ij and references therein, one should start not from 
space-time but from a symmetry algebra. Therefore in nonrelativistic quantum me- 
chanics we should start from the Galilei algebra and consider its irreducible represen- 
tations (IRs). For simplicity we again consider a one dimensional case. Let P^ = P 
be one of representation operators in an IR of the Galilei algebra. We can implement 
this IR in a Hilbert space of functions x{p) such that \x{p)\'^dp < oo and P is 
the operator of multiplication by p, i.e. Pxip) = VXip)- Then a question arises how 
the operator of the x coordinate should be defined. In contrast to the momentum 
operator, the coordinate one is not defined by the representation and so it should 
be defined from additional assumptions. Probably a future quantum theory of mea- 
surements will make it possible to construct operators of physical quantities from the 
rules how these quantities should be measured. However, at present we can construct 
necessary operators only from rather intuitive considerations. 

By analogy with the above discussion, one can say that semiclassical wave 
functions should be of the form xip) = exp{—ixop/h)a{p) where the amplitude a{p) 
has a sharp maximum near p = po & [pi,P2] such that a{p) is not small only when 
p G [pi, P2] • Then Ap is of the order of p2 —pi and the condition that the momentum is 
semiclassical is Ap ^ |po|. Since ihdx{p)/dp = Xoxip) + ihexp{—ixop/h)da{p)/dp, we 
see that x{p) be approximately the eigenfunction of ihd/dp with the eigenvalue 
Xq if |xoa(p)| ^ h\da{p) / dp\. Since \da{p)/dp\ is of the order of |a(p)/Ap|, we have 
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a condition |a;oAj9| ^ h. Therefore if the coordinate operator is X = ihd/dp, the 
uncertainty of coordinate Ax is of the order of h/Ap, \xq\ ^ Ax and this imphes 
that the coordinate defined in such a way is also semiclassical. We can also note that 
\xQAp\/2TTh is approximately the number of oscillations which the exponent makes on 
the segment [pi,P2] and therefore the number of oscillations should be much greater 
than unity. It is also clear that semiclassical approximation cannot be valid if Ap 
is very small, but on the other hand, Ap cannot be very large since it should be 
much less than pq. By analogy with the above discussion, the requirement that the 
operator ihd/dp becomes the coordinate in classical limit does not define the operator 
uniquely. In nonrelativistic quantum mechanics it is assumed that the coordinate is 
a well defined physical quantity even on quantum level and that ihd/dp is the most 
pertinent choice of the operator of this quantity. 

The above results can be directly generalized to the three-dimensional 
case. For example, if the coordinate wave function is chosen in the form 

V^(r) = ^ + ^Por] (3) 

then the momentum wave function is 

X(p) = y exp(--pr)^(r)^^-^ = -^,exp{ ^(p-po)ro] (4) 

It is easy to verify that 

ll^ir = 1 \^{^)fd\ = 1, llxir = / \xip)\'d'p = 1, (5) 

the uncertainty of each component of the coordinate operator is a/ ^2 and the un- 
certainty of each component of the momentum operator is h/{a\/2). Hence Eqs. ([3]) 
and (HI) describe a state which is semiclassical in a maximal possible extent. 

A well-known fact of quantum theory is that there is no operator having 
the meaning of the time operator. Hence a problem arises how time should be under- 
stood in quantum theory and this problem is discussed in a wide literature (see e.g. 
Ref. p]). It is usually assumed that time is a classical parameter such that the de- 
pendence of the wave function on time is defined by the Hamiltonian according to the 
Schroedinger equation. As noted in Sec. IH in some cases even such an interpretation 
of time might be problematic but in typical situations this assumption can be treated 
as a good approximation. In nonrelativistic quantum mechanics the Hamiltonian of 
the particle with the mass m is H = p^/2m and hence, as follows from Eq. (jlj), in 
the model discussed above the dependence of the momentum wave function on t is 
given by 

/ ,x Q^^^ r (p - Po)^a^ i. . ipH 

^(P' ^) = T^'''^^ h^^ - P<^)^° - 
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It is easy to verify that for this state the mean value of the operator p and the 
uncertainty of each momentum component are the same as for the state x(p); 
those quantities do not change with time. 

Consider now the dependence of the coordinate wave function on t. This 
dependence can be calculated by using Eq. ([6]) and the fact that 



/i dp 
expi-pr)x{p,t) ^^^^^3/2 (7) 



The result of a direct calculation is 

where = p^/mis the semiclassical velocity. This result shows that the semiclassical 
wave packet is moving along the classical trajectory r(t) = ro + vot. At the same 
time, it is now obvious that the uncertainty of each coordinate depends on time as 
Axj(t) = Axj(0)(l + /i^tVm2a^)i/2 (j = 1^2, 3) where Axj(O) = a/v^, i.e. the width 
of the wave packet in coordinate representation is increasing. This fact, known as the 
wave-packet spreading (WPS), is described in many textbooks and papers (see e.g. 
the textbooks [3J and references therein). It shows that if a state was semiclassical in 
the maximal extent at t = 0, it will not have this property at t > and the accuracy 
of semiclassical approximation will decrease with the increase of t. The characteristic 
time of spreading can be defined as = ma^/h. For macroscopic bodies this is 
an extremely large quantity and hence in macroscopic physics the effect of the wave- 
packet spreading can be neglected. In the formal limit % — > cxd, t^, becomes infinite, i.e. 
spreading does not take place. This shows that WPS is a pure quantum phenomenon. 

As shown in Ref. 0] titled " Nonspreading wave packets", for a one- 
dimensional wave function in the form of an Airy function, spreading does not take 
place and the maximum of the quantity |'?/'(x)p propagates with constant accelera- 
tion even in the absence of external forces. Those properties of Airy packets have 
been observed in optical experiments |5]. However, since such a wave function is not 
normalizable, we believe that the term "wave packet" in the given situation might be 
misleading since the mean values and uncertainties of the coordinate and momentum 
cannot be calculated in a standard way. Such a wave function can be constructed only 
in a limited region of space. As explained in Ref. this wave function describes 
not a particle but rather families of particle orbits. As shown in Ref. |4], one can 
construct a normalized state which is a superposition of Airy functions with Gaus- 
sian coefficients and "eventually the spreading due to the Gaussian cutoff takes over" . 
This is an additional argument that the effect of WPS is an inevitable consequence 
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of quantum theory. 

3 Wave packets for relativistic particles in momen- 
tum space 

The usual approach to Poincare symmetry on quantum level follows. Since classical 
Minkowski space is invariant under the action of the Poincare group, in quantum 
theory the angular momentum operators M'^'" and the four-momentum operators 
{fi, u = 0,1, 2, 3; M^'^ = —M'^^^) should satisfy the commutation relations of the 
Poincare group Lie algebra 

[M'"', MP""] = -i{7]PPM'''' + r]^''MPP - r]^''M''P - 7]''pM^'') (9) 

where rj^" is the diagonal metric tensor such that 7f^ = — r^^^ = —77^^ = —7]^^ = 1 and 
for simplicity Eqs. (Q are written in units h = c = 1. This approach is in the spirit 
of the well-known Klein's Erlangen program in mathematics. 

However, as we argue in Ref. [1] and references therein, quantum theory 
should not be based on classical space-time background and the approach should be 
the opposite. Each system is described by a set of independent operators. By defi- 
nition, the rules how these operators commute with each other define the symmetry 
algebra. In particular, by definition, Poincare symmetry on quantum level means that 
the operators commute according to Eq. (j9]). 

The next step in our construction is the definition of elementary particle. 
Although theory of elementary particles exists for a rather long period of time, there 
is no commonly accepted definition of elementary particle in this theory. In Ref. 
[T] and references cited therein we argue that, in the spirit of Wigner's approach 
to Poincare symmetry [6], a general definition, not depending on the choice of the 
classical background and on whether we consider a local or nonlocal theory, is that 
a particle is elementary if the set of its wave functions is the space of an irreducible 
representation (IR) of the symmetry algebra in the given theory. 

There exists a wide literature describing how IRs of the Poincare algebra 
can be constructed. In particular, an IR can be implemented in a space of functions 
x(p) such that the momentum operator P is the operator of multiplication by p. For 
particles with spin, the functions x(p) also depend on the spin projections but we 
will not write this dependence explicitly. 

As follows from Eqs. (Q, the operator I2 = — P^, where E = P°, is the 
Casimir operator of the second order, i.e. it is a bilinear combination of representation 
operators commuting with all the operators of the algebra. As follows from the well- 
known Schur lemma, all states belonging to an IR are the eigenvectors of I2 with the 
same eigenvalue which is denoted as m^. Then the energy operator is the operator 
of multiplication by ±e(p) where e(p) = (m^ -|- p^)^/^. The choice of the energy sign 
is only a matter of convention but not a matter of principle. Indeed, the energy can 
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be measured only if the momentum p is measured and then it is only a matter of 
convention what sign of the square root should be chosen. However, it is important 
that the sign should be the same for all particles. For example, if we consider a 
system of two particles with the same values of and the opposite momenta pi and 
P2 such that pi + P2 = 0, we cannot define the energies of the particles as e(pi) and 
— e(p2), respectively, since in that case the total four- momentum of the two-particle 
system will be zero what contradicts experiment. 

The notation I2 = rn^ is justified by the fact that for all known particles 
this quantity is greater or equal than zero. Then the mass m is defined as the square 
root of m? and the sign of m is only a matter of convention. The usual convention 
is that m > 0. However, from mathematical point of view, IRs with /2 < are not 
prohibited. If the velocity operator v is defined as v = P/i? then for known particles 
|v| < 1, i.e. |v| < c in standard units. However, for IRs with /2 < we will have 
that |v| > c and, at least from the point of view of mathematical construction of 
IRs, this case is not prohibited. The hypothetical particles with such properties are 
called tachyons and their possible existence is widely discussed in the literature. If 
the tachyon mass m is also defined as the square root of iv? then this quantity will be 
imaginary. However, this does not mean than the corresponding IRs are unphysical 
since all the operators of the Poincare group Lie algebra can depend only on m^. 

The usual choice of the Hilbert space H for an IR is such that the func- 
tions belonging to H are quadratically integrable over the Lorentz invariant measure 
(i^p/e(p). However, it is always possible to perform a unitary transformation such 
that the Hilbert space will be implemented as a space of functions quadratically in- 
tegrable over d^p. For reasons which will be clear below it is convenient for us to use 
such a normalization. 

We conclude that in relativistic quantum theory the four-momentum op- 
erators are well defined and have a clear physical meaning. In particular, it is possible 
to construct IRs describing the photon and other elementary particles. 

Consider first a construction of the wave packet for a particle with nonzero 
mass. A possible way of the construction follows. We first consider the particle in 
its rest system, i.e. in the reference frame where the mean value of the particle mo- 
mentum is zero. The wave function Xo(p) iii this case can be taken as in Eq. (jl]) 
with po = 0. As noted in the preceding section, such a state cannot be semiclas- 
sical. However, it is possible to obtain a semiclassical state by applying a Lorentz 
transformation to Xo(p)- As shown in a wide literature, in standard quantum theory 
(based on complex numbers) any IR representation of the algebra ^ by Hermitiam 
operators can be extended to an unitary IR of the Poincare group. One can show 
(see e.g. Eq. (2.4) in Ref. [7]) that for a spinless particle the unitary representation 
operator U{g) corresponding to a Lorentz transformation g can be defined as 

U{g)xM = [^]'^\o{p') (10) 
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where p' is the momentum obtained from p by the Lorentz transformation g ^. li g 
is the Lorentz boost along the z axis with the velocity v then 



P± = P±, Pz= (I _ ,,2)1/2 (11) 



where we use the subscript ± to denote projections of vectors onto the xy plane. 

As follows from this expression, exp(— p'^a^/2/l^) as a function of p has 
the maximum at p_L = 0, Pz = Pzo = v[{m'^ + p^^) /(I — v"^)]^^"^ and near the maximum 

a^p'^ 1 
— ^) - exp{-^[a^pl + b^Pz-Pzof]} 

where b = a(l — t>^)^/^ what represents the effect of the Lorentz contraction. If 
m ^ h/a (in units where c = 1) then m ^ |p^| and p^o ~ mv/{l — f^)^/^. In this 
case the transformed state is semiclassical and the mean value of the momentum is 
exactly the classical (i.e. nonquantum) value of the momentum of a particle with 
mass m moving along the z axis with the velocity v. However, in the opposite case 
when m <^ h/a the transformed state is not semiclassical since the uncertainty of pz 
is of the same order as the mean value of Pz- 

If the photon mass is exactly zero then the photon cannot have the rest 
state. However, even if the photon mass is not exactly zero, it is so small that 
the relation m <^ h/a is certainly satisfied for any realistic value of a. Hence a 
semiclassical state for the photon or a particle with a very small mass cannot be 
obtained by applying the Lorentz transformation to Xo(p) and considering the case 
when V is very close to unity. In this case we will describe a semiclassical state by a 
wave function which is a generalization of the function (jl]): 



a6^/^ Pio^ (Pz — Pofb"^ 



I 



Here we assume that the vector po is directed along the z axis and its z component is 
Pq. In the general case the parameters a and b defining the momentum distributions 
in the transverse and longitudinal directions, respectively, can be different. In that 
case the uncertainty of each transverse component of momentum is h/{a\/2) while 
the uncertainty of the z component of momentum is h/{b\/2). In view of the above 
discussion one might think that, as a consequence of the Lorentz contraction, the 
parameter b should be very small. However, the above discussion shows that the 
notion of the Lorentz contraction has a physical meaning only if m ^ h/a while 
for the photon the opposite relation takes place. We will see below that in typical 
situations the quantity b is large and much greater than a. 



4 The photon in semiclassical approximation 

Consider now a problem whether an elementary particle can be described in terms 
of space-time characteristics. By analogy with the discussion in Sec. [2], we can 
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state that in relativistic quantum theory there is no operator corresponding to time 
and hence at best t can be only a good parameter defining the evolution of the wave 
function according to the Schroedinger equation with the relativistic energy operator. 
A problem also arises whether it is possible to define a physical coordinate operator 
for an elementary particle. For example, in the representation where H is the space of 
functions quadratically integrated over d^p one might define the position operator by 
analogy with nonrelativistic quantum mechanics i.e. as ihd/dp. This is the Newton- 
Wigner position operator proposed in Ref. [8]. Then the coordinate wave function 
ip{r) is again defined by Eq. (jS]) and a question arises whether this operator has all 
the required properties of the physical coordinate operator. 

Let us first make a few remarks about the terminology of quantum theory. 
The terms "wave function" and "particle-wave duality" have arisen at the beginning 
of quantum era in efforts to explain quantum behavior in terms of classical waves but 
now it is clear that no such explanation exists. The notion of wave is purely classical; 
it has a physical meaning only as a way of describing systems of many particles by 
their average characteristics. In particular, such notions as frequency and wave length 
can be applied only to classical waves, i.e. to systems consisting of many particles 
such that space-time characteristics of those systems are measured on classical level. 
If a particle state vector contains exp[i{px — Et)/h] then by analogy with the theory of 
classical waves one might say that the particle is a wave with the frequency u = E/h 
and the (de Broglie) wave length A = 2'Kh/p. However, such defined quantities 
u and A are not real frequencies and wave lengths measured e.g. in spectroscopic 
experiments where only characteristics of many-particle systems are measured. In 
quantum theory the photon and other particles can be characterized by their energies, 
momenta and other quantities for which there exist well defined operators. Those 
quantities might be measured in collisions of those particles with other particles. The 
term "wave function" might be misleading since in quantum theory it defines not 
amplitudes of waves but only amplitudes of probabilities. So, although in our opinion 
the term "state vector" is more pertinent than "wave function" we will use the latter 
in accordance with the usual terminology. 

In classical theory the notion of field, as well as that of wave, is used for 
describing systems of many particles by their average characteristics. For example, 
the electromagnetic field consists of many photons. In classical theory each photon 
is not described individually but the field as a whole is described by the quantities 
E(r, t) and B(r, t) which can be measured (in principle) by using macroscopic test 
bodies. In particular, the notions of electric and magnetic fields of a single photon 
have no physical meaning. 

It has been well-known since the 1930s that, when quantum mechanics 
is combined with relativity, there is no operator satisfying all the properties of the 
spatial position operator. In other words, the coordinates cannot be exactly mea- 
sured even in situations when exact measurements are allowed by the non-relativistic 
uncertainty principle. For example, in the introductory section of the well-known 
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textbook [9] the following arguments are given in favor of this statement. Suppose 
that we measure the coordinates of an electron with the mass m. When the uncer- 
tainty of the coordinates is of the order of h/mc, the uncertainty of momenta is of 
the order of mc, the uncertainty of energy is of the order of mc^ and hence creation 
of electron-positron pairs is allowed. As a consequence, it is not possible to locahze 
the electron with the accuracy better than its Compton wave length h/mc. Hence, 
for a particle with a nonzero mass the exact measurement is possible only either in 
the non-relativistic limit (when c — )■ oo) or classical limit (when ^ — > 0). If m = 
is possible, the problem becomes even more complicated since the photon can create 
other photons with lesser energies. However, those arguments do not exclude a pos- 
sibility that the Newton- Wigner position operator can be meaningful in semiclassical 
approximation. 

In standard textbooks on quantum electrodynamics (see e.g. Ref. jlOj ) 
it is stated that in this theory there is no way to define a photon wave function 
in coordinate representation and the arguments are as follows. The electric and 
magnetic fields of the photon in coordinate representation are proportional to Fourier 
transforms of |p|^''^x(p), rather than x(p)- As a consequence, the quantities E(r) and 
B(r) are defined not by ip{r) but by integrals of ip{r) over some region. However, this 
argument also does not exclude the possibility that ip{r) can have a physical meaning 
in semiclassical approximation since, as noted above, the notions of the electric and 
magnetic fields of the single photon do not have a physical meaning. 

One of the arguments that the Newton- Wigner position operator does not 
have all the required properties follows. If at t = the function ip{r) has a finite 
carrier (i.e. ip{r) ^ only if r belongs to a vicinity of some vector Tq) and the 
evolution of ip{r,t) is governed by the Schroedinger equation with the relativistic 
energy operator then it is easy to show that at any t > the carrier of '?/'(r, t) will 
belong to the whole three-dimensional space. Then at any t > the particle can 
be detected at any point of the space and this contradicts the requirement that no 
information should be transmitted with the speed greater than c. However, one might 
say that the requirement that no signal can be transmitted with the speed greater 
than c has been obtained in Special Relativity which is a classical (i.e. nonquantum) 
theory operating only with classical space-time coordinates. As noted above, from the 
point of view of quantum theory the existence of tachyons is not prohibited. Note also 
that when two electrically charged particles exchange by a virtual photon, a typical 
situation is that the four-momentum of the photon is spacelike, i.e. the photon is the 
tachyon. 

A rather striking example is a photon emitted in the famous 21cm transi- 
tion line between the hyperfine energy levels of the hydrogen atom. The phrase that 
the lifetime of this transition is of the order of r = 10^ years implies that the width 
of the level is of the order of h/r, i.e. experimentally the uncertainty of the photon 
energy is h/r. Hence the uncertainty of the photon momentum is %/{ct) and with 
the above definition of the coordinate operators the uncertainty of the longitudinal 
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coordinate is cr, i.e. of the order of 10'' light years. Then there is a nonzero proba- 
bihty that immediately after its creation at point A the photon can be detected at 
point B such that the distance between A and B is 10^ light years. 

A problem arises how this phenomenon should be interpreted. On one 
hand, one might say that in view of the above discussion it is not clear whether or 
not the requirement that no information should be transmitted with the speed greater 
than c should be a must in relativistic quantum theory. On the other hand, as pointed 
out to me by Alik Makarov, we can know about the photon creation only if the photon 
is detected and when it was detected at point B at the moment of time t = to, this 
does not mean that the photon travelled from A to B with the speed greater than 
c since the time of creation has an uncertainty of the order of 10'' years. After the 
detection, the uncertainties of the coordinates have the order of the dimensions of 
the detector. Note also that in this situation a description of the system (atom + 
electric field) by the wave function (e.g. in the Fock space) depending on a continuous 
parameter t has no physical meaning (since roughly speaking the quantum of time 
in this process is of the order of 10'' years). If we accept this explanation then we 
should acknowledge that in some situations a description of evolution by a continuous 
classical parameter t is not physical. This is in the spirit of the Heisenberg S-matrix 
program that in quantum theory one can describe only transitions of states from the 
infinite past when t — )■ —oo to the distant future when t — )■ +oo. 

The above discussion shows that on quantum level the physical meaning of 
the coordinate is not clear but at least there are reasons to think that in some cases 
the transverse component of the Newton- Wigner position operator has a physical 
meaning in semiclassical approximation (see also Sec. |5]). If we also assume that 
in some situations time is a good approximate parameter describing the evolution 
according to the Schroedinger equation with the relativistic Hamiltonian then the 
dependence of the momentum wave function (fT2l) on t is given by 

X(p, t) = exp{--pct)x{p, 0) (13) 

where p = |p| and we assume that the particle is ultrarelativistic, i.e. p^ m. 

In view of the above discussion, the function ip{r,t) can be again defined 
by Eq. where now x(p, t) is defined by Eq. ( fT3l) . If the variable pz in the integrand 
is replaced by po + Pz then as follows from Eqs. ( 17|12fl3p 



-Yl(P^+Por + plY^'}d'p (14) 

We now take into account the fact that in semiclassical approximation the quantity po 
should be much greater than the uncertainties of the momentum in the longitudinal 
and transversal directions, i.e. po ^ pz and po ^ \p±\- Hence with a good accuracy 
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we can expand the square root in the integrand in powers of |p|/po- Taking into 
account the hnear and quadratic terms in the square root we get 

[{Pz + PoY + pi]'/' ^Po+Pz + pi/2po (15) 

Then the integral over d^p can be calculated as the product of integrals over (Pp± and 
dpz and the calculation is analogous to that in Eq. (j8]). The result of the calculation 
is 

^(r,t) = [n'/^ab'/'il + '-^)]-^exp[Upor - Pod)] 

(r± - ro±)^(l - ^) (^z-zo-cty 
^2 9.9 KT^ (16) 

This result shows that the wave packet describing an ultrarelativistic par- 
ticle (including a photon) is moving along the classical trajectory z{t) = zq + ct, 
in the longitudinal direction there is no spreading while in the transversal direction 
spreading is characterized by the function 



The characteristic time of spreading can be defined as t^: = p^a? /he. The fact that 
— 7- oo in the formal limit % ^ Q shows that in relativistic case WPS also is a pure 
quantum phenomenon (see the end of Sec. [2]) . From the formal point of view the result 
for is the same as in nonrelativistic theory but m should be replaced hy E/(? where 
E is the energy of the ultrarelativistic particle. However, in the ultrarelativistic case 
spreading takes place only in the direction perpendicular to the particle momentum. 
If t ^ the transversal width of the packet is a{t) = hct/poa. Hence the speed of 
spreading in the perpendicular direction is v^: = hc/poa. 



5 Geometrical optics 

The relation between quantum and classical electrodynamics is well-known and is de- 
scribed in textbooks (see e.g. Ref. [10]). As already noted, classical electromagnetic 
field consists of many photons and in classical electrodynamics the photons are not 
described individually. Instead, classical electromagnetic field is described by field 
strengths which represent average characteristics of a large set of photons. For con- 
structing the field strengths one can use the photon wave functions x(P)^) or ip{r,t) 
where E is replaced by hu and p is replaced by hk. In this connection it is interesting 
to note that since a; is a classical quantity used for describing a classical electromag- 
netic field, the photon is a pure quantum particle since its energy disappears in the 
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formal limit h ^ 0. Even this fact shows that the photon cannot be treated as a 
classical particle and the effect of WPS for the photon cannot be neglected. 

With the above replacements the functions x "wiH not contain any 

dependence on h (note that the normalization factor h'^^"^ in x(k, t) will disappear 
since the normalization integral for is now over (i^k, not d^p). The quantities 

u and k are now treated, respectively, as the frequency and the wave vector of the 
classical electromagnetic field and the functions x(k, t) and '?/'(r, t) are interpreted not 
such that they describe probabilities for a single photon but such that they describe 
classical electromagnetic field and E(r, t) and B(r, t) can be constructed from these 
functions as described in textbooks on quantum electrodynamics (see e.g. Ref. [ID]). 

As noted in the preceding section, some authors (see e.g. Ref. [10]) state 
that the function ^/'(r, t) cannot be interpreted as the photon wave function in coor- 
dinate representation since for each value of r, E(r, t) and B(r, t) depend not only 
on ip{r,t) but on the values of the function ip in some vicinity of r. This vicinity 
has dimensions of the order of the wave length. Hence, for example for visible light, 
where the wave length is of the order of hundreds of nanometers, the quantity r can 
be treated as a good coordinate in semiclassical approximation. Another argument 
in favor of this statement is that in classical electrodynamics the quantities E(r, t) 
and B(r, t) for the free field should satisfy the wave equation d'^'Ei/c^dt^ = AE and 
analogously for B(r, t). Hence if E(r,t) and B(r, t) are constructed from ip{r,t) as 
described in textbooks (see e.g. Ref. [ID]), they will satisfy the wave equation since, 
as follows from Eqs. fl7|12|13p . 4'{T,t) also satisfies this equation. 

The approximation of geometrical optics can be formulated in full analogy 
with semiclassical approximation in quantum theory. This approximation implies 
that if ko and tq are the average values of the wave vector and the spatial radius 
vector for a wave packet describing the electromagnetic wave then the uncertainties 
Ak and Ar, which are the average values of |k — ko| and |r — ro|, respectively, should 
satisfy the requirements Ak ^ |ko| and Ar ^ |ro|. Analogously, in full analogy 
with the derivation of Eq. ([2]), one can show that for each j = 1, 2, 3 the uncertain- 
ties of the corresponding projections of the vectors k and r satisfy the requirement 
AkjAvj > 1/2 (see e.g. Ref. [H]). In particular, an electromagnetic wave satisfies 
the approximation of geometrical optics in the greatest possible extent if AkAr is of 
the order of unity. 

We conclude that in the language of classical waves the parameters of 
spreading can be characterized by the function a{t) (see Eq. f[T7|) ) and the quantities 
and v^: (see the end of the preceding section) such that in terms of the wave length 
A = 2iTc/uJo 

In Ref. [12] the problem of WPS for classical electromagnetic waves has 
been discussed in the Fresnel approximation (i.e. in the approximation of geometrical 
optics) for a two-dimensional wave packet. Equation (25) of Ref. [12] is a special case 
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of Eq. ( !T5|) and the author of Ref. [12] shows that, in his model the wave packet 
spreads in the direction perpendicular to the group velocity of the packet. As noted 
at the end of the preceding section, in the ultrarelativistic case the function a{t) is 
given by the same expression as in the nonrelativistic case but m is replaced hj E/c^. 
Hence if the results of the preceding section are reformulated in terms of classical 
waves then m should be replaced by huo/c^ and this fact has been pointed out in 
Ref. Il2j. 

The quantity A^n = b/X shows how many oscillations the oscillating expo- 
nent in Eq. ( IT6ll makes in the region where the wave function or the amplitude of the 
classical wave is significantly different from zero. As noted in Sec. Ej for the validity 
of semiclassical approximation this quantity should be very large. In nonrelativistic 
quantum mechanics a and b are of the same order and hence the same can be said 
about the quantity N = a/X. As noted above, in the case of the photon we don't 
know the relation between a and b. In terms of the quantity N we can rewrite the 
expressions for t^, and in Eq. f[T5]l as 

U = 27iNX ^* = ^ (19) 

where T is the period of the classical wave. Hence the accuracy of semiclassical ap- 
proximation (or the approximation of geometrical optics in classical electrodynamics) 
increases with the increase of A^. 



6 Discussion 

In plain language, the problem discussed in this paper can be formulated as follows. 
Is the fact that we can see distant stars and even planets compatible with the another 
known fact that the wave function of the photon which has managed to survive its 
long journey to the Earth was the subject of WPS? 

We discussed WPS assuming that photons coming to the Earth do not 
interact with each other and hence WPS for each photon can be discussed indepen- 
dently. Is this assumption legitimate? For example, as shown in standard textbooks 
on quantum optics (see e.g. Ref. [13] and references therein)) quantum states de- 
scribing the laser emission are strongly coherent meaning that they are quantum 
superpositions of states with different numbers of photons for each value of momen- 
tum and polarization. In this case the approximation of independent photons is not 
legitimate. However, laser emission can be created only at very special conditions 
when energy levels are inverted, the emission is amplified in the laser cavity etc. At 
the same time, the main part of the radiation emitted by stars is understood such 
that it can be approximately described as the blackbody radiation and in addition a 
part of the radiation consists of photons emitted from different atomic energy levels. 
In that case the emission of photons is spontaneous rather than induced and one 
might think that the photons can be treated independently. Several authors (see e.g. 
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Ref. [13] and references therein) discussed a possibility that at some conditions the 
inverted population and amplification of radiation in stellar atmospheres might occur 
and so a part of the radiation can be induced. This problem is now under investi- 
gation. Hence we adopt a standard assumption that a main part of the radiation is 
spontaneous. 

The next question is whether the interaction of photons in the light emit- 
ted by stars is important or not. As explained in standard textbooks on quantum 
electrodynamics (see e.g. Ref. [IQj), the photon-photon interaction can go only via 
intermediate creation of virtual electron-positron or quar-antiquark pairs. If u is the 
photon frequency, m is the mass of the charged particle in the intermediate state and 
e is the electric charge of this particle than in the case when Tiu ^ mc^ the total 
cross section of the photon-photon interaction is [TO] 



a = ^^ — {—Y{^f (20) 
birm'^ 90^ he mc^ 

For photons of visible light the ratio huj/{mc^) is very small and a is an extremely 
small quantity. At present the effect of the direct photon-photon interaction has not 
been detected, and experiments with strong laser fields were only able to determine 
the upper limit of the cross section [15] . 

Our conclusion is that with a good accuracy the effect of WPS can be 
considered for each photon independently. Hence a question arises what can be said 
about the characteristics of photons coming to the Earth from distance objects. Typ- 
ical conclusions based on numerous experiments with light coming to the Earth from 
the Sun are as follows. We know that with a good accuracy this light can be de- 
scribed in the framework of geometrical optics i.e. one can approximately treat the 
light as a collection of particles moving along classical trajectories. Since we know 
that the photons came from the Sun then with a good accuracy we know the direction 
of their momenta and since we know the distribution of wave lengths then (in the 
approximation described in Sec. E]) we know the distribution of photon energies. 

The next question is what we know about the width of the coordinate 
photon wave functions in the direction perpendicular to the photon momentum in 
the approximation when the coordinate operators are defined as in Sec. |H Suppose 
that a wide beam of light falls on a screen which is perpendicular to the direction 
of light. Suppose that the total area of the screen is S but the surface contains 
slits with the total area Si. We are interested in the question of what part of the 
light will pass the screen. One might think that the obvious answer is that the part 
equals Si/S. This answer follows from the picture that the light consists of many 
photons moving along geometrical trajectories and hence only the Si/ S part of the 
photons will pass the surface. Numerous experiments show that deviations from the 
above answer begin to manifest in interference experiments where dimensions of slits 
and distances between them have the order of tens or hundreds of microns or even 
less. Hence one can conclude that the width of the photon wave functions cannot 
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be of the order of say centimeters or meters since in that case deviations from the 
Si/S law would be visible if the slits and the distances between them would have the 
corresponding dimensions but this does not happen. 

Consider, for example, the Lyman transition 2P — )■ IS* in the hydrogen 
atom on the Sun. In this case the energy of the photon is E = 10.2eV, its wave 
length is A = 121. 6nm, the lifetime is r = 1.6 ■ lO^^s and the period of the wave is 
T ^ 4 - 10~^^s. Since the lifetime is rather small, one might think that a description of 
the process by a continuous parameter t is a good approximation (see the discussion 
in Sec. H]). Hence the phrase that the lifetime is r can be interpreted such that the 
uncertainty of the energy is ti/r, the uncertainty of the longitudinal momentum is 
h/cT and b is of the order of cr ^ 0.48m or greater. In view of the above discussion, 
the estimation a ^ b seems to be very favorable since one might expect that the 
value of a is much less than 0.48m. With this estimation N = a/X ^ 4 ■ 10^. So the 
value of N is rather large and in view of Eq. f lT^ one might think that the effect of 
spreading is not important. However, this is not the case since, as follows from Eq. 
( IT9|) . ~ 0.04s. Since the distance between the Sun and the Earth is approximately 
t = 8 light minutes and this time is much greater than t^,, the width of the wave 
packet when it arrives to the Earth is v^,t ~ 5760m. It is obvious from the above 
discussion that such a value of the width is unrealistically large. On the other hand, 
if we assume that the initial value of a is of the order of several wave lengths then 
the value of N is much less and the width of the wave packet coming to the Earth is 
much greater. 

Consider now a photon which was created in the same reaction but on 
Sirius which is the brightest star on our sky. Since the distance to Sirius is 8.6 light 
years, an analogous estimation shows that even in the favorable scenario the width 
of the wave packet coming to the Earth from Sirius will be approximately equal to 
3 ■ lO^km but in less favorable situations the width will be much greater. 

As already noted, a standard understanding of light coming to the Earth 
from the Sun is such that the major part of the light comes not from transitions 
between atomic levels but from processes which can be approximately described as 
a black body radiation. In that case the spectrum of the radiation is approximately 
continuous and we cannot estimate the quantity a as above. However, even if we take 
for a a very favorable value of the same order as above we obviously will come to the 
same conclusion that the width of the wave packet will be unreasonably high. 

One might say that the expression for in Eqs. flTK]) and f lTI?]) resembles 
a well-known phenomenon of diffraction: if a wave encounters an obstacle having a 
dimension d it begins to diverge and the angle of divergence is of the order of X/d. 
However, the phenomenon of WPS implies that the width of the wave packet in the 
transversal direction is growing even when the packet propagates in empty space. This 
phenomenon takes place only for wave function in coordinate representation while the 
distribution of momenta remains unchanged. As a consequence, even when at some 
moment of time the packet was maximally semiclassical (i.e. for each component of 
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the coordinate and momentum the product of their uncertainties is of the order of 
h), this property is not conserved with time since uncertainties of coordinates in the 
transversal direction become greater. 

The problem of WPS in the ultrarelativistic case has been discussed in a 
wide literature. As already noted in the preceding section, in Ref. [12] the effect of 
WPS has been discussed in the Fresnel approximation for a two-dimensional model 
and the author shows that in the direction perpendicular to the group velocity of 
the wave spreading is important. He considers WPS in the framework of classical 
electrodynamics. We believe that considering this effect from quantum point of view 
is even simpler since the photon wave function satisfies the relativistic Schroedinger 
equation which is linear in d/dt. As noted in Sec. [5l this function also satisfies 
the wave equation but it is simpler to consider an equation linear in d/dt than that 
quadratic in 8/ dt. However, in classical theory there is no such an object as the photon 
wave function and hence one has to solve either a system of Maxwell equations or the 
wave equation. Probably the approach of Ref. [12] can be generalized such that the 
results of the present paper can be recovered. However, in Ref. [12] the effect of WPS 
is discussed in view of diffraction and interference experiments with electromagnetic 
waves and electrons but not in view of photons coming to the Earth from distant 
objects. There is also a number of works where the authors consider WPS in view 
of propagation of classical waves in a medium such that dissipation is important (see 
e.g. Ref. [16]). In a recent paper [17] the effect of WPS has been discussed in view of 
a possible existence of superluminal neutrinos. The authors of this work consider only 
the dynamics of the wave packet in the longitudinal direction in the framework of the 
Dirac equation. They conclude that wave packets describing ultrarelativistic fermions 
do not experience WPS in this direction. However, the authors do not consider WPS 
in the perpendicular direction. 

We believe that in view of many works discussing the effect of WPS it 
is rather strange that (to the best of our knowledge) a natural question whether or 
not WPS is important for photons which can travel to the Earth even for billions 
of years has not been raised in the literature. As already noted, probably the main 
reason is a belief that the effect of WPS is not important for such photons since we do 
not observe "fuzziness" in the location of astronomical objects. Also many physicists 
believe that it is a good approximation to treat photons coming to the Earth from 
distant objects as classical particles moving along classical trajectories, for examples 
as bullets. However, the phenomenon of WPS takes place for bullets too. In this 
sense the difference between bullets and photons is only quantitative: since bullets 
have large masses and sizes, the characteristic time of spreading for them is extremely 
large (see the end of Sec. [2]) while for photons this time is much smaller. 

In view of the above discussion a problem arises why we can see separate 
stars at all. Indeed, if the width of a wave packet is growing as discussed in the above 
examples then we would see not separate stars but only a continuous background of 
light coming to us from many stars. To the best of our knowledge, this problem state- 
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ment has not been discussed in the hterature. Below we discuss possible approaches 
for solving the problem. 

One of the possibilities might be such that the interaction of light with the 
interstellar medium cannot be neglected. On quantum level a process of propagation 
of photons in the medium is rather complicated because several mechanisms of prop- 
agation should be taken into account. For example, a possible process is such that a 
photon can be absorbed by an atom and reemitted. This process makes it clear why 
the speed of light in the medium is less than c: because the atom which absorbed the 
photon is in an excited state for some time before reemitting the photon. However, 
this process is also important from the following point of view: even if the coordinate 
photon wave function had a large width before absorption, as a consequence of the 
phenomenon known as the collapse of the wave function, the wave function of the 
emitted photon will have in general much smaller dimensions since after detection 
the width is defined only by parameters of the corresponding detector. If the photon 
encounters many atoms on its way, this process does not allow the photon wave func- 
tion to spread significantly. Analogous remarks can be made about other processes, 
for example about rescattering of photons on large groups of atoms, rescattering on 
elementary particles if they are present in the medium etc. However, such processes 
have been discussed in Theory B and, as noted in Sec. [H they probably result in 
more blurring that is seen. 

The interaction of light with the interstellar medium might also be im- 
portant in view of hypotheses that the density of the interstellar medium is much 
greater than usually believed. Among the most popular scenarios are dark energy, 
dark matter etc. As shown in our papers (see e.g. Ref. [1] and references therein), the 
phenomenon of the cosmological acceleration can be easily and naturally explained 
from first principles of quantum theory without involving dark energy, empty space- 
background and other artificial notions. However, the other scenarios seem to be more 
realistic and one might expect that they will be intensively investigated. A rather 
hypothetical possibility is that the propagation of light in the interstellar medium 
has something in common with the induced emission when a photon induces emission 
of other photons in practically the same direction. In other words, the interstellar 
medium amplifies the emission as a laser. This possibility seems to be not realistic 
since it is not clear why the energy levels in interstellar matter might be inverted. 

A fully different possibility might be as follows. As discussed in Sees. 
I2MI while the choice of the momentum operator is unambiguous, the choice of the 
coordinate operator is not. We described the coordinate photon wave function by 
using the Newton- Wigner position operator which in turn is a generalization of the 
standard position operator in nonrelativistic quantum mechanics. The validity of 
this operator has been proved by many experiments in atomic and nuclear physics. 
However, as argued in Ref. [18], for describing the motion of macroscopic bodies 
on astronomical distances standard position operator should be modified. As shown 
in this paper (see also Ref. [T]), it is possible to modify the position operator such 
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that the Newton law of gravity and the precession of Mercury's perihehon will be 
recovered. A problem arises whether the effect of WPS will take place if the motion 
of photons is described by the modified rather than standard position operator. This 
problem is very important and deserves a serious investigation. 

We conclude that WPS of photons coming to the Earth from astronomical 
objects is the fundamental unsolved problem of physics. An interesting historical 
analogy with this situation follows. In 1935 Hubble and Tolman comparing Theory A 
and Theory B wrote in Ref. [19] that they "... both incline to the opinion, however, 
that if the red-shift is not due to recessional motion, its explanation will probably 
involve some quite new physical principles". Probably, an analogous conclusion is 
valid in the present situation. 
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